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ABSTRACT 

We present an analysis of X-ray and ultra-violet data of the dwarf nova VW 
Hyi that were obtained with XMM-Newton during the quiescent state. The X-ray 
spectrum indicates the presence of an optically thin plasma in the boundary layer that 
cools as it settles onto the white dwarf. The plasma has a continuous temperature 
distribution that is well described by a power-law or a cooling flow model with a 
maximum temperature of 6-8 keV. We estimate from the X-ray spectrum a boundary 
layer luminosity of 8 x 10'^° erg s~^, which is only 20 per cent of the disk luminosity. 
The rate of accretion onto the white dwarf is 5 x 1O~^^M0 yr"-'^, about half of the 
rate in the disk. From the high-resolution X-ray spectra, we estimate that the X-ray 
emitting part of the boundary layer is rotating with a velocity of 540kms~^, which 
is close the rotation velocity of the white dwarf but significantly smaller than the 
Keplerian velocity. We detect a 60-s quasi-periodic oscillation of the X-ray flux that 
is likely due to the rotation of the boundary layer. The X-ray and the ultra-violet flux 
show strong variability on a time scale of ^--^1500 s. We find that the variability in the 
two bands is correlated and that the X-ray fluctuations are delayed by ~100 s. The 
correlation indicates that the variable ultra-violet flux is emitted near the transition 
region between the disk and the boundary layer and that accretion rate fluctuations 
in this region are propagated to the X-ray emitting part of the boundary layer within 
~100 s. An orbital modulation of the X-ray flux suggests that the inner accretion disk 
is tilted with respect to the orbital plane. The elemental abundances in the boundary 
layer are close to their solar values. 

Key words: accretion, accretion discs - binaries: close - stars: dwarf novae - stars: 
individual: VW Hyi - novae, cataclysmic variables - X-rays: stars 



1 INTRODUCTION 

Dwarf novae are powerful X-ray sources with luminosities 
of lO'^^-lO"^'^ ergs~^. The origin of these X-rays is thought 
to be the boundary layer between the white dwarf and the 
inner edge of the accretion disk. There the disk material, ini- 
tially moving at a Keplerian velocity, dissipates its kinetic 
energy while being decelerated to the rotation velocity of 
the white dwarf. The structure of the boundary layer is still 
an unsolved problem in our understanding of equatorial ac- 
cretion onto white dwarfs. In the standard picture of dwarf 
novae, half of the accretion energy is radiated away by the 
disk, primarily at optical and ultra-violet (UV) wavelengths. 
Unless the white dwarf is rotating rapidly, the other half of 
the energy should be liberated in the boundary layer and 
emitted as X-ray and extreme ultra-violet (EUV) photons. 
However, observations have shown that for many dwarf no- 



vae the boundary layer luminosity is much lower than the 
disk luminosity. In VW Hyi the X-ray flux during quies- 
cence was found to b e a factor of 4 lowe r than expected in 
the standard picture iBelloni et al.llT99lh . 

At a distance of 65 pc JWarnerlligSTT) . VW Hyi is one 
of the nearest and brightest cataclysmic variables. It has 
been studied extensively at optical, UV, and X-ray wave- 
lengths. VW Hyi is a dwarf nova of type SU UMa and un- 
dergoes normal outbursts every ~30 d and superoutbursts 
every ~180 d. It has an orbital period of 107 min and lies 
in a direction with an extr emely low hydrogen column den- 
sity Nh = 6x 10"cm-2 dPolidan et alJll990tl . The X-ray 
spectrum in quiescence is consistent with emission from an 
optically thin, multi-temperature plasma with temperatures 
imto -12keV dWheatlev et alJll99dl . lvan der Woerd et alJ 
il987fl detected coherent soft X-ray oscillations during su- 
peroutburst with a period of 14 s and suggested that they 
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originate from the difFerentially rotat inR outer layers of the 
white dwarf or the boundary layer. ISion et all i200ll) in- 
ferred that the white dwarf's rotation velocity projected 
along the line of sight is 400-500 kms~^, corresponding to 
a rotation period of ~80 s. Fo r our calculations, w e assume 
an orbital inclination i = 60° ( iRltter fc Kolblll99?j). a white 
dwarf mass M^d = 0.63 Mg JSchoembs fc VoidllQsil) . and 
a white dwarf radius R^d = 8.3 x 10^ cm ( derived from the 
mass-radius relationship in lHamada fc S jpeter 19( jjl. 

In this paper we present our analysis of XMM-Newton 
data obtained while VW Hyi was in a quiescent state. With 
the high sensitivity of the X-ray detectors and their broad 
coverage of the X-ray band, we are able to determine the 
temperature distribution of the plasma in the boundary 
layer. We search the X-ray and UV light curves for periodic 
and quasi-periodic oscillations and characterize the proper- 
ties of the low- frequency noise (flickering). The capability of 
XMM-Newton for simultaneous X-ray and UV observations 
enables us to study the correlation between the variability of 
the boundary layer and the disk luminosity. From the high- 
resolution RGS spectrum, we infer the rotational velocity of 
the X-ray emitting plasma in the boundary layer. We derive 
abundances for elements with strong X-ray emission lines. 



2 OBSERVATIONS AND DATA REDUCTION 

VW Hyi was observed with XMM-Newton l|,Tansen et alJ 
I2OOII) on 2001 October 19. At the time of the observation, 
VW Hyi was in a quiescent state 22 d after a normal out- 
burst and 23 d before a superoutburst (derived from light 
curves provided by the AAVSO). Froni the EPIC instru- 
ments (|Turner et al.ll200ll : IStriider et alll200j) we obtained 
18.7 ks (MOS) and 16.1 ks (PN) of continuous X-ray data. 
The EPIC cameras were operating in small window mode 
(MOS) and in full frame mode (PN). For MOS-1 and PN 
the medium blocking filters and for MOS -2 the thin blocking 
filter were selected. The Optical Monitor iMason et all200lr) 
performed 5 exposures with a total duration of 16.1 ks using 
the UVWl filter (240-340 nm) and one exposure of 1.2 ks 
with the B filter. The Optical Monitor (OM) was operated 
in fast mode with a sample time of 0.5 s. We ob tained 19.3 ks 
of data with the RGS ( den Herder et ai]|200lll . Background 
levels were low during the entire observation. 

From the EPIC MOS and PN data, we extracted source 
photons using a circular aperture with a radius of 40 arc- 
sec. The count rates given in this paper have been corrected 
for the 88 per cent enclosed energy fraction of this aperture 
relative to a 60 arcsec aperture (see XMM-Newton Users' 
Handbook) . The background rate, which was estimated from 
larger off-center regions, contributes less than 1 per cent to 
the count rate in the source aperture. We included in our 
analysis good photon events (FLAG=0) with patterns 0- 
12 for MOS and 0-4 for PN. To create X-ray light curves, 
we applied a barycentric correction to the photon arrival 
times and combined the events from all three EPIC instru- 
ments. From the OM fast mode data, we extracted source 
photons using a circular aperture of 5 arcsec radius (81 per 
cent enclosed energy fraction). Sky background contributed 
~2 per cent to the count rate in the source aperture. The 
high UV brightness of VW Hyi caused significant coinci- 
dence losses at a level of 10-20 per cent. We corrected for 



these coincidence losses using the method described in the 
XMM-Newton Users ' Handbook. 



3 TIMING ANALYSIS 

3.1 X-ray and UV light curves 

Fig. shows background subtracted X-ray and UV light 
curves covering three orbital cycles of the binary (orbital 
period 107.0 min). As orbital phase we selected BJD(TT) 
2452201.748, which is the time of peak optical bright- 
ness as predicted by the ephemeris in Ivan Amerongen et alJ 
(1987). The ephemeris has an accumulated uncertainty of 
250 s or 0.04 in phase. With this choice of orbital phase 
0, inferior conjunction of the secondary occurs at ~0.27 
JSchoembs fc Vogdll98ll) . 

The X-ray light curves in Figs. and |5| show variability 
on two time scales. Clearly visible is an orbital modulation 
of ~30 per cent that is somewhat anti-correlated with the 
UV flux. The hardness ratio in Fig. |21 indicates that the X- 
ray spectrum is softer during the fainter part of the orbital 
cycle. The low- frequency flickering common in dwarf novae 
is seen as strong flaring with amplitudes up to ~40 per cent 
and a recurrence time scale of ~1500 s (see Section l3.3|l . 
Despite the strong variability of the X-ray flux on this time 
scale, the hardness ratio in Fig. remains fairly constant. 
The average count rate in the 0.2-12 keV energy range is 
4.0 s-\ 

The UV light curve in Figs. and |5| show a strong or- 
bital modulation in the form of a broad hump. This or- 
bital hump is thought to be due to emission from the bright 
spot, the hot region where the accretion stream f rom the 
secon dary impacts the outer edge of the disk (e.g. IWarnerl 
I1995I) . On time scales shorter than the orbital period, the 
UV flux is less variable than the X-ray flux (see Section f3.3l 
for a statistical analysis). This is to be expected since the 
UV flux is dominated by emission from the inner accretion 
disk and the white dwarf, while the X-rays mostly originate 
from the boundary layer. 

In the UVWl fllter, a count rate of 1 corresponds 
to a flux of 4.5 X 10"^*^ erg cm" ^ s"^A"^ at 290 nm. Ac- 
cordingly, the UV flux from VW Hyi was varying between 
2.5-3.8 X 10"^'* erg cm~'^ s~^ with an average value of 
2.9 X 10"^'' ergcm"^ s"^A"^ or 8.1 mjy. From the short 
B-band exposure (not shown), we estimate that the B mag- 
nitude was varying between 14.4 and 13.9. Both the UV 
and the B-band brightness ranges agree with past measure- 
ments for VW Hyi during quiescence dPringle et al.lir987l : 
Ivan Amerongen et al.llT98^ . 

3.2 Periodic and quasi-periodic oscillations 

Periodograms of the X-ray and UV light curves are shown in 
Fig. El The spectra l power is given in Leahy normalization 
jLeahv et alJll983l) . which is suitable for data from pho- 
ton counting detectors. In this normalization, the power in 
each spectral bin has a ^'^-distribution with 2 degrees of 
freedom (mean 2 and variance 4), provided that the source 
is non-variable. The dotted lines in Fig. |21 are the 90 per 
cent confidence levels for the detection of strictly periodic 
signals. For frequencies above 0.003 Hz, we do not find any 
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Figure 1. X-ray and UV light curves of VW Hyi obtained with XMM-Newton on 2001 October 19 (see Section 1,3. li . Also shown is the 
ratio of the count rates in the 1-12 keV and 0.2-1 keV energy ranges. VW Hyi has an orbital period of 107.0 min. 
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Figure 2. X-ray/UV light curves and X-ray hardness ratio folded 
on the orbital period of 107.0 min. 



evidence for periodic signals. Although the detection thresh- 
old is exceeded below 0.003 Hz, the large spectral power is 
not due to periodic signals but rather the stochastic process 
of flickering common in dwarf novae (see Section l8.3|l . For 
the frequency range 0.005-1 Hz, we derived upper limits of 
3.3 per cent and 0.9 per cent, respectively, for the amplitude 



of any periodic signal in the X-ray and UV data (confidence 
level 90 per cent). 

To search for quasi-periodic signals, we divided the light 
curves into shorter time intervals and averaged the peri- 
odograms for all intervals. The insets in Fig.|3show the av- 
erages of 10 such periodograms. For the X-ray data, there is 
a clear excess of spectral power between 0.01-0.02 Hz. Since 
the power is spread over many frequency bins, the excess 
is not due to a strictly periodic signal but rather indicates 
the presence of a quasi-periodic oscillation (QPO). By sum- 
ming the spectral power in the frequency range 0.011-0.021 
Hz, we find that the QPO is statistically significant at a 
6(7 level. To estimate the parameters of the QPO, we fitted 
to the periodogram a spectral density model consisting of 
a Gaussian plus a constant term for the photon counting 
noise using a maximum likelihood method. We obtained a 
center frequency of 0.016 Hz (62 s period), a FWHM fre- 
quency spread of 0.006 Hz (53-77 s period range), and an 
oscillation amplitude (RMS) of 6.2 per cent. To determine 
the quality of the fit, we simulated spectra with the power 
density function obtained from the fit by randomizing am- 
plitudes and phases according to the appropriate probabil- 
ity distributions. We find that 96 per cent of the simulated 
power spectra yield a worse fit than the XMM-Newton data. 

To further characterize the QPO, we assume that the 
peak in the power spectrum is caused by a series of random 
pulses in the light curve that have a Gaussian profile with 
a sinusoidal modulation. From the above fit results, we can 
then derive some average properties of these pulses. We ob- 
tain a pulse duration of 110 s (FWHM of the Gaussian pro- 
file) and Ap/tp — 0.004s~'^, where Ap is the pulse amplitude 
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Figure 3. Periodograms of the X-ray and UV light curves shown with logarithmic (left) and linear (right) scales. The spectral power 
is given in Leahy normalization for which each bin has an average value of 2 provided that the source is non-variable and that the 
observational noise is only due to photon counting statistics (see Section lij. 21 . The dotted lines indicate the 90 per cent confidence levels 
for the detection of strictly periodic signals and the dashed line marks the orbital frequency Fprb- The solid curves show the best fit 
of a power-law spectral density function. Averaged periodograms obtained by splitting the light curves into 10 equal time intervals are 
shown in the insets. 



in counts-s"^ and tp the average pulse separation. Ap and tp 
cannot be determined independently, but the non-detection 
of individual pulses in the light curve requires Ap < 1 s~^ 
and therefore tp < 250 s. 

The QPO is detected both for hard (>1 keV) and soft 
(<1 keV) X-rays with no apparent difference in strength. 
However, the QPO is not detected in the UV light curve, 
which places an upper limit of 0.7 per cent on the RMS 
amplitude in the 0.011-0.021 Hz range. Our periodogram 
analysis did not reveal any other QPOs in the frequency 
range 0.02-1 Hz. This places an upper limit of 2.4 per cent 
for the X-ray data and 0.7 per cent for the UV data on the 
RMS amplitude of any QPO with a frequency spread of 0.01 
Hz. 

3.3 Flickering and X-ray/UV correlation 

At low frequencies, the power spectrum is dominated by 
the flickering typical for dwarf novae. This low-frequency 
variability is commonly modelled with a power-law spectral 
densi ty function S{f) oc /"'' (see e.g. iTimmer fc Koeni3 
Il995^ . The power-law index f3 is indicative of the underlying 
process. For example, random walk noise has an index — 2. 
We fitted the periodograms with this power-law model plus a 
constant component for the photon counting noise. To avoid 
a bias from the orbital variations at a frequency of 0.000156 
Hz, we restricted the fit to frequencies above 0.0002 Hz. 
The results of the fits are shown in Fig. |H] (solid curves). We 
obtained for the power-law index /3 a value of 1.90±0.19 for 
the X-ray data and 2.02±0.14 for the UV data. Both indices 
are consistent with a random walk noise process. To quantify 
the strength of the variability, we quote the RMS amplitude 
of the power-law noise above 0.0002 Hz. We obtained an 



RMS of 18 per cent for the X-ray data and 7 per cent for the 
UV data. To determine the quality of the fits, we simulated 
power spectra as described in Section f3.2l We found that in 
X-rays 45 per cent and in the UV 93 per cent of the simulated 
spectra yield a worse fit than the XMM-Newton data. We 
investigated the dependence of the power-law parameters on 
the X-ray energy but could not find statistically significant 
differences between hard and soft X-rays. 

At frequencies below ~0.0005 Hz, the X-ray spectral 
power appears to be lower than predicted by the model fit, 
indicating that the power-law does not continue toward low 
frequencies. We determined a 99 per cent probability that 
the low spectral power in the range 0.0002-0.0005 Hz is 
not a coincidence and inconsistent with the power-law fit. 
If frequencies below 0.0005 Hz are excluded from our fit, the 
power-law index increases to 2.2 ± 0.2, which is still consis- 
tent with P — 2. The concentration of spectral power around 
0.0007 Hz, which gives rise to the strong variability on a 
~1500-s time scale (Fig. 0, may indicate the presence of a 
QPO. However, the duration of the XMM-Newton observa- 
tion is insufficient to distinguish a QPO from a power-law 
spectral density function with a changing index 

To investigate the connection between the X-ray and 
the UV variability, we calculated the cross-correlation be- 
tween the two light curves. We removed the obvious cor- 
relation on the orbital time scale by first applying a high- 
pass filter with a cut-off frequency at 0.0004 Hz. The result- 
ing cross-correlation coefficient as a function of time lag is 
shown in Fig. 2] A strong peak near zero time lag suggests a 
significant correlation between the X-ray and the UV vari- 
ability. The correlation coefficient is normalized such that it 
has a value of 1 for maximally correlated signals. The high- 
frequency photon-counting noise introduces a bias that can 
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Figure 4. Cross-correlation between X-ray and UV light curves 
calculated after the variability on time scales longer than 2500 s 
had been removed (see Section lii.ijI . The offset of the central peak 
with respect to indicates that the X-ray variations lag behind 
those in the UV by ~100 s. 



reduce the correlation coefficient. If we correct for this bias, 
the peak correlation coefficient increases to 0.5. The central 
peak is offset with respect to zero time lag, indicating that 
the X-ray variations are delayed relative to those in the UV 
by ~100 s. The X-ray/UV correlation disappears when the 
cut-off frequency of the high-pass filter is increased to 0.001 
Hz, demonstrating that the correlation is indeed due to the 
low-frequency flickering on the 1500-s time scale. We also 
investigated the correlation between the hard and soft X- 
ray variability. We found a strong correlation between the 
X-ray fluxes above and below 1 keV, but the time lag was 
consistent with zero to within a few seconds. 



4 SPECTRAL ANALYSIS 

The medium resolution EPIC spectrum and the high reso- 
lution RGS spectrum are shown in Figs. |S] and |S] Both spec- 
tra exhibit strong emission lines, indicating the presence of 
a hot, optically thin plasma. The prominent lines identifled 
by our model flts have been labelled. Spectral fit ting was 
perfo rmed with the XSPEC package version 11.2 (|A rnaudl 
|l99y). As spectral models we used multi-temperature ver- 
sions of the MEKAL model, which describes an o ptically thin, 
isothe rmal plasma based on c alculations by iMewe et alJ 
1^8^ and lLiedahl et al] il995l) . For optimal fit results, we 
binned the EPIC and RGS spectra at one third of the respec- 
tive FWHM detector resolution. In addition, a minimum of 
20 counts per bin were required for the EPIC spectra so that 
X^-statistics could be used. To retain the high resolution of 
the RGS spectra, we did not further rebin the RGS data but 
instead accounted for the low number of counts per bin by 
using C-statistics (ICashlll979l) . 



4.1 EPIC spectrum 

The X-ray spectra of dwarf novae in quiescence are com- 
monly fit with one or two single-temperature plasma models. 
We attempted to fit these models to the EPIC spectra but 
found that they are inconsistent with the data. However, 
a good fit could be achieved with a combination of three 



single-temperature mekal models. For the temperatures of 
the three components we found 0.5, 1.5, and 5keV. Adding a 
fourth component to the model did not improve the fit. It is 
evident from this result that the X-ray emitting plasma has a 
continuous temperature distribution, at least over the range 
0.5-5keV. The multi-temperature nature can be expected for 
accreting gas that cools as it settles onto the white dwarf. 
To account for the spread in temperature, we attempted 
to fit several MEKAL-based continuous-temperature models 
provided by XSPEC. 

The best results were achieved with the CEMEKL model, 
which fit the data slightly better than the three-temperature 
model but with a smaller number of parameters. As the 
simplest of the MEKAL-based multi-temperature models, CE- 
MEKL assumes that the Emission Measure (EM) has a 
power-law temperature-dependence 

where Tmax is the maximum temperature of the plasma and 
DEMiTmax) the Differential Emission Measure at Tmax 



DEM{Tmax 



dEM 
d logT 



(2) 



The integration in the CEMEKL code is performed logarith- 
mically (note the log T above), so that a con stant EM{T) 
corresponds to a = 1 JPone fc OsborndllQQ'il) . The results 
of our fits are summarized in Table Q Rows 1-3 show the 
parameters obtained when the spectrum from each of the 
three EPIC cameras was fit individually. As can be seen in 
the null hypothesis probabilities, the CEMEKL model is in 
very good agreement with the data. The three EPIC spec- 
tra yield somewhat different values for the parameters a 
and Tmax- We attribute this to small discrepancies in the 
cross-calibration of the EPIC detectors. Because the param- 
eters a and Tmax are strongly anti-correlated, the differences 
between the three fits are not as significant as they might 
appear. Row 4 in Table shows the result for a simulta- 
neous fit of the three EPIC spectra. The null hypothesis 
probability for the combined fit is significantly lower than 
for the individual fits, which is a further indication for the 
cross-calibration discrepancies. Assuming isotropic emission 
and a distance of 65 pc, we obtain from the fit a Differ- 
ential Emission Measure DEM{Tmax) = 4.3 x lO'^^cm"^ 
at kTmax ~ 6.0 keV. Integrating the power-law model from 
T = to Tmax yields a total Emission Measure EM — 
1.4 X lO^'^cm^''. The neutral hydrogen colu mn density Nh 
towa rd VW Hyi is fairly low at 6 x lO'^'^cm"^ jPolidan et alJ 
I1990I) . Accordingly, adding photoelectric absorption to our 
model did not improve the fit. We derived an upper limit 
for Nh of 2.2 x 10^** cm'^ (at 90 % c.L). 

In the CEMEKL model, the elemental abundance ratios 
are fixed at solar values and only the total abundance is var- 
ied. A significant improvement of the fit could be achieved 
with the CEVMKL model (solid curve in Fig. |K| and row 6 
in Table 0. In this model, the abundances for the 13 most 
common elements from C to Ni can be varied independently. 
Table |5| shows the best-fit values and the 90 per cent confi- 
dence ranges of the abundances obtained with the CEVMKL 
model. Good constraints are found for the elements O, Mg, 
Si, S, and Fe, all of which have abundances near the solar 
values. The abundances are mainly derived from the K-a 
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Figure 5. X-ray spectra obtained with the EPIC PN and MOS-1 cameras. Prominent K-a emission lines have been labeled. A magnified 
view of the Fe K-a line is shown in the inset. For clarity, the MOS-1 data points have been shifted down by a factor of 2. The solid 
curves show the best fit with the CEVMKL model (see Section l4.1l and Table [Tl . 



Table 1. Results of our spectral fits of various optically thin, multi-temperature plasma models to the EPIC data (see Section l4.1i . The 
CEMEKL and CEVMKL models assume a power-law temperature distribution of the Emission M easure. MKCFLOW is a mod el for a cooling 
flow. Except for the CEVMKL model, the abundance ratios were fixed at the solar values given in lAnders fc Grevessd ^^S^)- For the latter 
three fits, the spectra from all three EPIC cameras were fit simultaneously. In the MKCFLOW model, the fitted low-temperature cut-off 
kTmin was 0.08 keV (minimum allowed by model). Uncertainties are given at a Icr level. 



Data sets 


XSPEC 


Power-law 


kTrnax 


Abundance 


Flux 


erg 


cm ^ s ^] 




(d.o.f.) 


Null hypothesis 




model 


index a 


[keV] 


[solar] 


0.2-12 keV 


bolometric 




probability 


PN 


CEMEKL 


1.24 ±0.04 


6.4 ±0.3 


1.03 ± 0.05 


6.4 X 10 


-12 


8.2 X 10" 


12 


212.1 (185) 


8.4 % 


MOS 1 


CEMEKL 


1.63 ±0.09 


5.3 ± 0.3 


0.99 ± 0.08 


6.1 X 10 


-12 


7.4 X 10~ 


12 


280.6 (169) 


25.7 % 


MOS 2 


CEMEKL 


1.38 ±0.05 


5.8 ± 0.3 


1.05 ± 0.08 


6.4 X 10" 


-12 


8.0 X 10- 


12 


174.0 (171) 


42.2 % 


PN-I-MOS 


CEMEKL 


1.35 ±0.03 


6.0 ±0.2 


1.02 ± 0.04 


6.3 X 10^ 


-12 


8.0 X 10" 


12 


605.3 (533) 


1.6 % 


PN-fMOS 


MKCFLOW 




7.8 ± 0.2 


0.97 ±0.04 


6.4 X 10 


-12 


7.8 X 10- 


12 


616.7 (533) 


0.7 % 


PN-fMOS 


CEVMKL 


1.28 ±0.04 


6.2 ± 0.2 


see Table 13 


6.3 X 10 


-12 


8.1 X 10- 


12 


557.3 (521) 


13.1 % 



emission line of each element. The most prominent of these 
lines are labeled in Fig. |^ Because of the lower spectral res- 
olution of the EPIC detectors, confusion of emission lines 
may bias the measurements for the less abundant elements. 
In particular, the estimates for Ne and Na are unreliable, 
since their K-a lines lie in the range 0.7-1.3 keV, which is 
dominated by Fe L-shell lines. 

A good fit, though slightly worse than with the CE- 
MEKL model, was achieved with the mkcflow model (row 
5 m Table m. MKCFLOW is based on the cooling flow model 
bv lMushotzkv fc SzvmkowiakI il988l) and assumes that the 
Emission Measure for each temperature is proportional to 
the inverse of the bolometric luminosity (i.e. proportional 
to the cooling time). The variable parameters in the model 
are the minimum and maximum temperatures of the flow 
{Tmin and Tmax), the elemental abundance (abundance ra- 
tios fixed at solar values), and the flux normalization. The 



good agreement of the X-ray spectrum with the mkcflow 
model indicates the presence of hot plasma in the boundary 
layer that cools as it settles onto the white dwarf. How- 
ever, the cooling flow model assumes constant pressure and 
spherical symmetry, which may not be valid in the strong 
gravitational field of the white dwarf and for the shear and 
turbulence in a rotating, belt-like boundary layer. Since the 
CEVMKL model provides the best fit, we will use it in our 
further analysis. 

The maximum plasma temperature of 6-8 keV in 
VW H yi is low compared to other dwarf novae. For in- 
stance, j^^^^^^I] (|200.?l ) found maximum temperatures 
of ^-^20 keV for U Gein and ~80 keV for SS Cyg, while 
iDone fc Osbornd lll997l) obtained ~20 keV for SS Cyg. 
From earlier ob s ervati ons of VW Hyi during qu iescence, 
iHartmann et all l)l99!j ) and lWheatlev et alJ (|l99(il) derived 
similarly low temperatures of 10 keV and 11 keV, respec- 
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Table 2. Elemental abundances derived from our spectral fits 
with the OEVMKL model (see Section s |4JJ and |4^ . Abundance s 
are given relative to the solar values in lAnders fc Grevessj flQSg^. 
Shown are the best-fit values and the 90 per cent confidence 
ranges. The more reliable abundance estimates have been high- 
lighted. 



Element 




EPIC 




RGS 




best fit 


90 % c.l. 


best fit 


90 % c.l. 


C 


1.5 


0.3-2.9 


0.82 


0.52-1.65 


N 


1.1 


0.0-2.4 


1.6 


0.9-2.9 


O 


0.95 


0.79-1.12 


0.80 


0.67-0.91 


Ne 


0.24 


0.0-0.55 


0.81 


0.29-1.26 


Na 


> 10 


> 3.4 


0.0 


0.0-10.4 


Mg 


1.27 


0.91-1.65 


0.9 


0.4-2.3 


Al 


4.5 


1.6-7.6 


0.0 


0.0-6.5 


Si 


1.30 


1.05-1.57 


2.8 


1.6-3.8 


S 


1.10 


0.77-1.45 


1.9 


0.9-2.9 


Ar 


0.0 


0.0-0.75 


4.8 


0.2-9.4 


Ca 


1.8 


0.6-2.9 


0.0 


0.0-3.0 


Fe 


1.10 


0.99-1.22 


0.68 


0.59-0.77 


Ni 


2.0 


1.1-3.0 


0.0 


0.0-1.1 



tively, using a cooling flow model. During these observations, 
however, the X-ray flux was about twice that observed with 
XMM-Newton (see below), which may explain our somewhat 
lower result of 7.8 keV. 

According to our model fits, the bolometric flux from 
the boundary layer is 8 x 10~^^ erg cm~^ s~^. Since most of 
the flux is emitted inside the XMM-Newton detector band, 
this estimate is fairly independent of the spectral model 
(see Table 0. For a distance of 65 pc, the correspond- 
ing boundary layer luminosity is Lu = 8 x 10"''' ergs^^. 
Here we assumed that, because of the high orbital incli- 
nation, half of the belt-like boundary layer is obscured 
by the white dwarf. We ignored reflection of X-rays by 
the white dwarf and a possible obscuration of the bound- 
ary layer by the disk. In comparison with previous ob- 
servations during quiescence, the X-ray flux observed with 
XMM-Newton was abo ut half of t hat found with Bepp oSAX 
jHartmann et alJll999ri. ROSAT llBelloni et al.lll99]J) . and 
EXOSAT jvan der Woerd fc Heis elll987h. and about twice 
that seen with Einstein llCord^v^^nTl9 8l]h . Despite the 
differences in X-ray flux, VW Hyi had a similar optical 
brightness during these observations. 

To further investigate the orbital modulation of the X- 
ray flux (Figs. and I^J, we extracted individual spectra 
for the bright and the faint parts of the hght curve (orbital 
phases 0.4-0.9 and 0.9-0.4, respectively). We find that the 
spectrum is softer during the faint phase and that the ratio 
of faint-phase to bright-phase spectrum is decreasing by ~10 
per cent from the lowest to the highest X-ray energies. The 
softening of the spectrum can also be seen in the hardness 
ratio in Fig. |5| This trend is opposite of what would be 
expected if the decline in brightness were due to absorption 
by some intervening gas. Possible causes for the orbital X- 
ray modulation are discussed in Section f5.4l 



4.2 RGS spectrum 

The RGS spectrum (Fig. ^ shows a number of prominent 
emission lines, in particular the Lyman-a lines of some com- 
mon elements and the Fe L-shell lines between 0.7 and 

1.3 keV. As for the EPIC spectra, we fitted the RGS spec- 
trum with the model for an optically thin plasma with a 
power- law temperature distribution (cevmkl). Because of 
the smaller energy range covered by the RGS, the model 
parameters a and Tmax were not well constrained. We there- 
fore fixed these parameters at the values obtained from the 
EPIC spectra (a — 1.28 and kTmax ~ 6.2 keV) and only 
fit the elemental abundances. The best-fit model is shown 
as a solid curve in Fig. |S| The elemental abundances (Table 
121 are in general agreement with those obtained from the 
EPIC spectra. The RGS spectrum tends to yield good es- 
timates for elements with prominent Lyman-a lines in the 
RGS detector band (C, O, Ne). For other elements (Mg, Si, 
S), the Lyman-Q lines are fairly weak and the more sen- 
sitive EPIC detectors provide better constraints. In Table 
121 the more reliable abundance estimates have been high- 
lighted. For most elements, our results are consistent with 
near solar abundances. We cannot confirm the l arge o ver- 
abundance of C and Si found by iHuang et alJ il996l) for 
a UV emitting ring rotating with the Keplerian velocity. 
However, the UV spectra on which these results are based 
were obtained 10 days after a superoutburst and the ring 
may ac tually be a post-o utburst equatorial belt on the white 
dwarf i Warner fc Woud t, ,2002) . The abundances in such a 
belt may be enhanced due to mixing with the outer layers 
of the white dwarf. In contrast, our estimates derived from 
the X-ray spectrum provide the abundances in the accreting 
gas. 

Good estimates from both the RGS and the EPIC spec- 
tra were obtained only for the abundances of oxygen and 
iron. While for oxygen the two values are consistent, there 
is a noticeable discrepancy for iron. In the EPIC spectrum, 
the iron abundance is mainly derived from the Fe xxv K- 
Q line at 6.7 keV, while in the RGS spectrum the L-shell 
lines of Fe xvil-xxiv determine the abundance. We interpret 
the discrepancy for iron as an indication that the power-law 
model slightly overestimates the emission measure for tem- 
peratures at which iron is in the ionization states Fe xvil- 
XXIV. The model appears to overpredict the strength of the 
Fe XX line at 12.8 A and the Fe xvil line at 15.0 A. The 
cause for this discre pancy may be a h i gh opa city for these 
lines. According to iBhatia fc Kastneil lll999ll . the Fe xvil 
line at 15.0 A weakens significantly relative to the other vis- 
ible Fe XVII lines for column densities >10^^ cm~^. For the 
accretion rate determined in Section l^).3l and for a boundary 
layer covering 10 per cent of the white dwarf, such a column 
density can be present if the iron in the cooling gas remains 
in the Fe xvil ionization state for ~5 s. 

The width of the O vill hue at 19.0 A is shghtly larger 
than the spectral resolution of the RGS. The likely cause 
is Doppler-broadening due to the rotation of the boundary 
layer around the white dwarf. Broadening due to the or- 
bital motion of the white dwarf itself is about one order of 
magnitude smaller and can be neglected. We fitted the line 
profile assuming that the X-ray emitting part of the bound- 
ary layer is a thin, equatorial belt near the white dwarf's 
surface and found a significantly non-zero rotation velocity 
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Figure 6. Flux-calibrated RGS spectrum showing the combined I''* and 2"'* order RGS-1 and RGS-2 data. Prominent emission lines 
have been labeled. The solid curve shows the best fit of the CEVMKL model (see Section l4.2l and Table [Tt. The data is binned at 0.072 A, 
the FWHM detector resolution. 



Vm sini — 520 km with a 90 per cent confidence range 
of 80-810 kms~^. Here i is the orbital inclination of the bi- 
nary. A slightly better constraint can be obtained by fitting 
the entire RGS spectrum as described earlier in this section 
but with the inclusion of the Doppler-broadening. We find 
a best-fit velocity Vu sini = 540 km s~^ with a 90 per cent 
confidence range of 280-810 kms~\ This resuH is m agree- 
ment with the rotation velocity of 750 km s~^ derived from 
the 60-s oscillation (see Section f5.1ll . 



5 DISCUSSION 

5.1 60-s X-ray oscillation and boundary layer 
rotation 

Coherent dwarf nova oscillations (DNOs) and quasi-periodic 
oscillations (QPOs) have been observed in many dwarf no- 
vae a t X-ray, ultra-violet and optical wavelength ( Warneil 
I1995I) . In VW Hyi the only X-ray oscillation previously 
found was a coherent 14 -s DNO during superoutburst 
ivan der Woerd et al.lll98^ . As its origin, the authors sug- 
gested the differentially rotating outer layers of the white 
dwarf or the boundary layer. At optical wavelength, DNOs 
have been observed in the period range 20-36 s and QPOs 
at the periods 23 s, 88 s, 253 s, an d 413 s ( s ee Ta ble 8.2 
in|Warncr 199^). Recently, Woudt fc Warneij (|20o3) found 
DNOs with a period evolving thoughout the outburst from 
14 s to >40 s and QPOs with evolving periods of hundreds 
of seconds. In all cases, the oscillations were seen during or 
shortly after an outburst. 

The 60-s X-ray QPO we found with XMM-Newton is 
the first oscillation observed in VW Hyi during quiescence. 
The QPO is detected with equal strength in the soft and 
hard X-ray bands (i.e. the spectral slope is similar to that of 
the overall X-ray spectrum), which indicates that the oscilla- 
tion is caused by small brightness variations in the boundary 
layer. Emission from a magnetic pole can be ruled out since 
the oscillation is not strictly periodic. The 60-s oscillation 
may be related to the rotation of the X-ray emitting gas in 
the boundary layer. Since VW Hyi is a high-inclination sys- 



tem {i = 60°), intermittently brightened parts of the bound- 
ary layer may be periodically obscured as they rotate around 
the white dwarf, thus creating a weakly coherent oscillation. 
The average live-time of the brightness fluctuations is ~110s 
(see Section r3.2L similar to the coherence time of the QPO. 

Assuming that the 60-s oscillation is due to rotation 
around the white dwarf, we can estimate the rotation veloc- 
ity of the X-ray emitting gas in the boundary layer. For an 
orbital inclination i = 60° and a radius equal to the white 
dwarf radius R^d = 8.3 x 10* cm, we find a rotation velocity 
Vbi sin i = 750kms~^ . This is slightly larger than the velocity 
of 540kms~^ we determined from the line-broadening (Sec- 
tion but still within the 90 per cent confidence range. 
The velocity derived from the QPO may be biased since 
it depends on the somewhat uncertain values of i and R^d- 
Our velocity measurements indicate that the boundary layer 
is rotating significantly slower than the inner accretion disk, 
which has a Keplerian velocity of ~3200 kms~^. However, 
the boundary layer velocity is close to the rotation velocity 
of the white dwarf V^^d sini = 400-500 km s~^ determined 
bv lSion et alj i200ll) from the width of UV-detected atmo- 
spheric lines. It appears that the accreting gas dissipates 
most of its rotational kinetic energy and decelerates to al- 
most the white dwarf's rotation velocity before it can radiate 
efficiently in X-rays. 

5.2 X-ray/UV correlation 

As discussed in Section 13.31 we found that the variability 
of the X-ray and UV fiux seen on a time scale of ~1500 s 
is correlated and that the X-ray fiux variations are delayed 
by ~100 s. In the past, correlation studies have been rare 
because of the technical difficulties of simultaneous X-ray 
and optical/UV observations. With the Optical Monitor 
on board XMM-Newton, the search for correlations is now 
greatly simplified. A similar correlati on as the one in VW 
Hyi was found bv lJensen et alJ il983l) in the nova- like vari- 
able TT Ari. There the X-ray variability was also delayed 
by ^1 min. 

The observed correlation is somewhat surprising since 



© 2003 RAS, MNRAS OOO.ITHT^ 



X-ray and UV observations of the dwarf nova VW Hyi in quiescence 9 



the X-ray and UV emission originate from distinct regions 
in the binary. While the source of the X-ray emission is the 
boundary layer very close to the white dwarf's su rface, the 
UV f lux is emitted by the inner accretion disk (e.g. lLa DousI 
Il98g^ and the white dwarf. Because of the relatively low 
temperatures, the disk or white dwarf do not contribute to 
the X-ray flux. Similarly, the contribution of the boundary 
layer to the UV flux, as extrapolated from the X-ray spec- 
trum, is only ~3 x 10~* of the total flux. 

Because the X-ray flux variations are trailing those in 
the UV, we can rule out reprocessing of X-rays in the ac- 
cretion disk as the cause of the correlation. Further investi- 
gations of the X-ray data did not revealed noticeable differ- 
ences between the spectrum during flares and the spectrum 
between flares. We therefore suspect that the X-ray variabil- 
ity is caused by fluctuations of the accretion rate onto the 
white dwarf. These fluctuations may be due to accretion rate 
variations in the disk, directly seen via the variable UV flux, 
that propagate inward to the boundary layer. However, the 
viscous time scale, i.e. the time scale for matter to diffuse 
inward through the d isk, is much longer than the observed 
delay of ~100 s (e.g. iFrank et al1ll992h . The variable UV 
flux must therefore be emitted from near the transition re- 
gion between the disk and the boundary layer. There the 
accreting gas is decelerated well below the Keplerian veloc- 
ity and is no longer supported by rotation against radial 
in-fall. Accretion rate fluctuations near the inner disk edge 
may be q uickly propagated to the inner par ts of the bound- 
ary laver. lMever fc Mever-Hofmeisterl ^ll994^ suggested that, 
during quiescence, accretion onto the white dwarf does not 
occur directly from the disk but via a corona that is formed 
by evaporation of the inner disk. In this model, the UV 
variability may be connected to variations of the evapora- 
tion rate that modulate the amount of gas in the corona 
and, with some time delay, the accretion rate onto the white 
dwarf. 

As a feasibility test, we estimate the minimum size of 
the region responsible for the variable UV flux. We assume 
that the region is emitting as a blackbody and that it is 
causing the 7 per cent variability of the flux in the UVWI 
filter (see Section FS.St . A further constraint on the blackbody 
can be derived from the non-de tection of UV em ission at 
wavelength shorter than 98 nm llLong et alJll996H . We find 
a maximum temperature of 3 eV and a minimum emitting 
area of 6 x 10^'' cm^. The corresponding width of an annular 
region with a radius equal to the white dwarf radius R^d is 
only 0.01 Rjnd, which is much smaller than the distance from 
the inner disk edge to the white dwarf (see below). Only a 
thin annular region near the inner disk edge is needed to 
produce the variable UV flux. 

If our interpretation of the X-ray/UV correlation is 
correct, the 100-s time delay can provide an estimate for 
the inner disk truncation radius. The radial in-fall velocity 
in the boundary la yer is given by K ~ a cJ/Vjc (see e.g. 
iNaravan fc Y^1994^ . where a is the viscosity parameter, 14 
the Keplerian velocity, and Cs = {kT/ finip)^^^ the isother- 
mal sound speed. With the assumption that the temperature 
in the boundary layer is roughly equal to the initial temper- 
ature of the cooling gas kTmax ~ 7keV fSection l4.1ll . and us- 
ing 14 = 3200 kms~^, a = 0.1, and a mean molecular weight 
H — 0.6, we find for the in-fall velocity Vr ^ 35 km s"'^. The 
100-s time delay then suggests that the inner disk edge is 



^OARwd above the white dwarf. This result depends on the 
rather uncertain viscosity parameter a and should only be 
considered a rough estimate. However, the inner truncation 
radius cannot be much larger than our estimate since this 
would be inconsistent with the observed orbital modulation 
of the X-ray flux (see Section [5.41 . 

5.3 X-ray luminosity and accretion rate 

In quiescence, the combined optical and UV flux from VW 
Hyi is 1.7 x 10~^° ergc m~'^s~^, h alf of which is probably 
due t o the accretion disk iPringle et al. 1987; Vorbunt et all 
[^S^). For a distance of 65 pc and an inclination of 60°, the 
corresponding disk luminosity is Ldisk = 4 x lO^'^ erg s~^. In 
Section [4. II we estimated for the boundary layer luminosity 
Lm = 8x 10'^° erg s~^. In agreement with earlier observations 
(e.g. iBelloni et al1ll99ll) , the ratio of boundary layer to disk 
luminosity Lu /Ldisk = 0.2 is well below the expected value 
of 1. The cause for this lack of boundary layer emission is still 
an open question in our understanding of dwarf novae. For 
VW Hyi we can rule out a rapidly rotating white dwarf as 
the cause since the rotation velocity of the boundary layer is 
much smaller than the Keplerian velocity (see Section f5. 11 . 

As discussed in Section 14.11 the X-ray spectrum is in 
good agreement with the model of a multi-temperature 
plasma having a cooling flow or power-law emission mea- 
sure distribution. This strongly suggests that the origin of 
the X-ray emission is cooling plasma in the boundary layer 
settling onto the white dwarf. The plasma has an initial tem- 
perature of 6-8 keV when it begins to radiate efficiently in 
X-rays. Our results are in rema rkable qualitative agr e ement 
with the model calculations bv lNaravan fc PophamI lll993h 
for dwarf novae at low accretion rates. The calculations pre- 
dict the presence of a narrow transition region between disk 
and boundary layer in which the accreting gas is decelerated 
well below the Keplerian velocity. At this inner disk edge, 
the gas begins to fall quickly toward the white dwarf and the 
density decreases by several orders of magnitude. Because 
of the low density, the gas can no longer cool efBciently and 
becomes invisible to us. The accreting material continues to 
dissipate its rotational kinetic energy while being heated to 
temperatures in excess of 10*K. Finally, when the in-falling 
gas is compressed as it piles up on the white dwarf, it begins 
to cool strongly via bremsstrahlung and becomes visible in 
X-rays. This compressed. X-ray emitting gas is localized in 
a thin region on the white dwarf and should be well approx- 
imated by a simple cooling flow model. In agreement with 
our measurements of the boundary layer rotation velocity 
(Section l5.11 . the model calculations predict that the gas in 
this thin region has been slowed down to almost the white 
dwarf rotation velocity. However, the model does not explain 
the low ratio Lm / Ldisk ■ 

As the in-falling gas is compressed and cools from an 
initial temperature Tmax to the much lower white dwarf tem- 
perature, it releases an enthalpy H = 5/2 x kTmax per parti- 
cle. We can therefore estimate the rate of accretion from the 
boundary layer onto the white dwarf Mti using the relation 

Lbl = — kTmax (3) 

(e.g. iFabianllTggi) . Here /i is the mean molecular weight 
(typically ~0.6) and rrip the proton mass. With an initial 
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temperature kTmax ~ 7keV and a boundary layer luminosity 
Lu = 8 X 10'^'^ ergs~^, we find an accretion rate Mu ~ 
5 X 10"^^Mo yr"\ 

For comparison, we estimate the accretion rate in the 
disk Mdisk from the disk luminosity L^isk- With the reason- 
able assumption that half of the accretion energy is dissi- 
pated in the disk, we can calculate Mdisk using the relation 



Ldisk — 



'2Rwd 



(4) 



For a white dwarf mass M^d = 0.63 and a radius 
R-wd = 8.3 X 10* cm, the disk accretion rate is Mdisk = 
12 X IQ-^^Mgyr-^ Our resuh seems to indicate that Mu ~ 
1/2 Mdisk, i.e. that half of the accreting gas is lost in a wind. 
A model for such an outflow in dwarf novae during qu i- 
escence is described in iMever fc Mever-HofmeisteJ il994l) . 
However, our estimate of Lu and therefore Mu may be too 
low if part of the boundary layer is obscured by the disk. The 
orbital modulation of the X-ray flux (Section 15.41 indicates 
that at least some obscuration is present. A geometrically 
thin disk can obscure at most half of the boundary layer 
so that the actual Mu could be up to twice our estimate. 
It is therefore still possible that Mu = Mdisk and that no 
outflow is present. Our estimate of Mdisk may be too low as 
well if the disk is truncated at a radius significantly larger 
than R-wd- In Section fS. 21 we found that the inner disk edge 
is probably ~0.4 R-^d above the white dwarf. This would 
imply that the actual Mdisk is ~40 per cent larger than 
our estimate. However, it seems unlikely that Mdisk can be 
sufficiently large so that the low ratio Lu/ Ldisk can be at- 
tributed to a low Mu/ Mdisk- In order for Mdisk to satisfy 
Mu/ Mdisk = Lu/ Ldisk = 0.2, the inner disk would have 
to be truncated at least ~1.0 R^d above the white dwarf 
(or more if we underestimated Mu)- However, such a large 
truncation radius is inconsistent with the observed orbital 
modulation of the X-ray flux (Section 15.41 . We conclude 
that the low boundary layer luminosity Lu cannot be solely 
attributed to a low accretion rate Mu- 

The key to understanding the missing boundary layer 
emission in VW Hyi may be the low initial temperature of 
the X-ray emitting gas. If the kinetic energy of the material 
at the inner disk edge, which is rotating at a Keplerian ve- 
locity of ~3200kms~^, is converted completely into thermal 
energy, the plasma would reach a temperature of ~21 keV. 
Yet when the in-falling gas first emits X-rays as it piles up on 
the white dwarf, its temperature is only 6-8 keV. The weak 
Fe XXVI line at 6.9keV (see Fig.lSJ clearly demonstrates that 
plasma at higher temperatures is either not present or not 
radiating efficiently in X-rays. It appears that the gas in the 
boundary layer looses ~2/3 of the available accretion energy 
before it begins to emit X-rays. 

It is possible that the gas inside the inner disk edge, 
where it is decelerated well below the Keplerian velocity, is 
cooling so rapidly that it cannot reach the maximum tem- 
perature of ~21keV. The gas would likely have temperatures 
higher than the inner disk but not high enough to emit X- 
rays that could be detected by XMM-Newton- However, no 
strong EUV emission from VW Hyi is observed, and, ow- 
ing to the low hydrogen column density, not much of the 
EUV flux can be hidden by interstellar absorption. From 
the non-dete ction of EUV em ission at wavelength shorter 
than 98 nm iLong et al.lll99^ and from the quiescent flux 



detected by EXOSAT dWheatlev et alJll996|) . we estimate 
that at most 20 per cent of the expected boundary layer 
lum inosity is emitted in the EUV. 

iLivio fc Pringld lll992fl found that a weak magnetic 
field (-10* G) of the white dwarf can truncate the in- 
ner disk, thus causing the UV delay observed during out- 
burst in some dwarf novae including VW Hyi. The den- 
sity inside the inner disk edge may be sufficiently low for 
the magnetic field to pentrate into the boundary layer and 
cause the plasma to emit cyclotron radiation. We estimate 
the cyclotron cooling time scale tcyd for thermalized, non- 
relati vistic electrons from the em itted cyclotron power (equ. 
6.7 in lRvbicki fc Lightma n 1979) and find t^y^ [s] = 4 x 10** x 
B [G]~^. For a magnetic field B = 10'* G, the cyclotron cool- 
ing time is ~4 s, comparable to the in-fall time of gas not 
supported by rotation. A significant fraction of the thermal 
energy could be lost to cyclotron radiation that would be 
emitted in the radio band. However, strong radio emissio n 
from dwarf novae is not observed (e.g. lBenz fc Guedell989ll . 
Also, we did not detect strictly periodic X-ray oscillations 
that would indicate omission from a magnetic pole. 

^ever & Movcr-Hofmeistcr ( 1994) suggested that dur- 
ing quiescence the inner accretion disk evaporates and forms 
a corona. The gas in the corona is partially accreted onto the 
white dwarf and partially lost in a wind. The X-ray spectrum 
is dominated by emission from the cooling plasma piling up 
on the white dwarf. In this model, some of the thermal en- 
ergy in the corona is lost to the wind or conducted back to 
the accretion disk, which may explain the low temperature 
of the plasma when it begins to emit X-rays. 



5.4 Orbital X-ray modulation 

Our investigation of the orbital X-ray modulation showed 
that the spectrum becomes softer when the fiux is decreas- 
ing (Section 14.11 . This spectral change is opposite of what 
would be expected if the decline in X-ray fiux were due to 
absorption by some intervening gas (e.g. an extended bright 
spot). We therefore suggest that the orbital X-ray modula- 
tion is caused by a geometric effect, i.e. that during the faint 
phase some hot regions of the boundary layer are blocked 
from view, while cooler regions remain visible. For this to 
be possible, the obscuring material has to be comparable in 
size and close to the boundary layer (e.g. the white dwarf or 
the inner accretion disk). We consider three possible expla- 
nations: 

(i) The boundary layer may not be axisymmetric. It is 
difficult to imagine how an asymmetry that is synchronized 
with the orbital motion can exist on the rapidly rotating 
white dwarf. One possibility is asymmetric accretion caused 
by a non-axisymmetric inner disk. Again, it is unclear how 
such an asymmetry in the disk could persist for several or- 
bital cycles. In some dwarf novae, stationary structures hav e 
been observed in the form of spiral waves fe.g. lSteegh!j200iri . 
albeit only during outbursts and in the outer parts of the 
disk. 

(ii) The boundary layer may be partially obscured by a 
non-axisymmetric inner disk. Since the transition from op- 
tically t hick disk to optically thin boundary layer is quite 
abrupt iNaravan fc Pophaml ll993'). the disk essentially acts 
like a circular aperture, blocking part of the X-ray emission. 
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The fraction of visible boundary layer should remain con- 
stant, yet it may vary with orbital phase if the inner edge of 
the disk is not circular. Again, it is unclear what the cause 
of such an asymmetry might be. 

(iii) The rotation axes of the white dwarf and the inner 
disk may not be aligned. It is thought that the boundary 
layer forms a belt-like structure on the white dwarf. Since 
the X-ray emitting part of the boundary layer is mostly co- 
rotating with the white dwarf, it is likely aligned with the 
equator. If the inner disk is tilted with respect to the white 
dwarf's equator, the fraction of boundary layer obscured by 
the disk may vary with phase, causing an orbital modulation 
of the X-ray flux. 

In the latter scenario, it is unlikely that the rotation 
axis of the white dwarf is tilted with respect to the orbital 
plane since any such m isalignment would disappear on a 
time scale of 100 years llKopaJll97dl . However, warping of 
the accretion disk may cau se a tilt of the inner disk edge. 
iMurrav fc Armitagel (^§3) found that tidal interactions are 
too weak to cause warping during the high state, but a tilt 
may b e pr e sent i n the quiescent disk of short-period dwarf 
novae. iLail l)l99!j^ suggested that warping of the inner disk 
may be caused by a weak magnetic field of the primary. The 
precession of a warped disk in VW Hyi could be revealed if 
future observations show a phase shift of the orbital X-ray 
modulation relative to the orbital hump seen in the UV. 



flux is well described by a power-law spectral density func- 
tion with an index consistent with a random walk noise pro- 
cess {P — 2). In X-rays the variability is caused by fluctua- 
tions of the accretion rate onto the white dwarf. We found a 
correlation between the X-ray and the UV variability with 
the X-ray fluctuations being delayed by ~100 s. This re- 
sult suggests that the variable UV flux is originating from 
near the transition region between disk and boundary layer 
and that accretion rate fluctuations in this region are prop- 
agated to the X-ray emitting part of the boundary layer 
within ~100 s. We estimated from the 100-s time delay that 
the accretion disk is truncated ~0.4 7?„d above the white 
dwarf. An orbital modulation of the X-ray flux suggests that 
the inner accretion disk is tilted relative to the white dwarf's 
equator and the orbital plane. This warping of the inner disk 
may be due to a weak magnetic field of the white dwarf. 
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6 CONCLUSIONS 

We have analyzed XMM-Newton data of the dwarf nova 
VW Hyi obtained during quiescence. The spectrum indi- 
cates that the X-ray emission is due to an optically thin 
plasma that cools as it settles onto the white dwarf. The 
cooling plasma has an initial temperature of 6-8 keV when 
it first emits X-rays. This temperature is well below the 
~21 keV expected if all the rotational kinetic energy of the 
material in the inner disk were instantly converted into heat. 
Apparently, the accreting gas looses ~2/3 of its thermal en- 
ergy before it piles up on the white dwarf and begins to 
radiate efficiently in X-rays. We estimated for the luminos- 
ity of the boundary layer 8 x 10^° erg s~^, which is only 20 
per cent of the accretion disk luminosity. The low boundary 
layer luminosity cannot be explained by a low mass transfer 
rate. We derived for the accretion rate onto the white dwarf 
5xlO"^^M0yr"^ and for that in the disk 12 x lO"^^M0yr"^ 
We determined elemental abundances for some common el- 
ements and found no large deviations from the solar values. 

From the Doppler-broadening of emission lines, we de- 
termined that the X-ray emitting part of the boundary 
layer is rotating with a velocity Vu sin i = 540 kms~^, 
which is similar to the white dwarf rotation velocity but 
much smaller than the Keplerian velocity of the inner disk. 
This result rules out a rapidly rotating white dwarf as the 
cause of the low boundary layer luminosity. We found a 60- 
s quasi-periodic X-ray oscillation with a coherence time of 
only ~100 s. The oscillation is likely due to the rotation 
of the boundary layer. The corresponding rotation velocity 
Vbi sin i = 750 km s~^ is consistent with that derived from 
the line-broadening. 

The low-frequency variability of the X-ray and the UV 
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